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Abstract 

Recurrence networks are complex networks, constructed from time series data, hav¬ 
ing several practical applications. Though their properties when constructed with 
the threshold value e chosen at or just above the percolation threshold of the net¬ 
work are quite well understood, what happens as the threshold increases beyond the 
usual operational window is still not clear from a complex network perspective. The 
present Letter is focused mainly on the network properties at intermediate-to-large 
values of the recurrence threshold, for which no systematic study has been performed 
so far. We argue, with numerical support, that recurrence networks constructed from 
chaotic attractors with e equal to the usual recurrence threshold or slightly above 
cannot, in general, show small-world property. However, if the threshold is further 
increased, the recurrence network topology initially changes to a small-world struc¬ 
ture and finally to that of a classical random graph as the threshold approaches the 
size of the strange attractor. 
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1 Introduction 


In the last two decades, complex network theory has emerged as a popular 
tool for analyzing complex and spatially extended systems mm- It has found 
applications in a wide range of helds including sociology [¥115] , communication 
P7] and biological sciences [8119] . The theory of complex networks initially 
started with random graphs (RG) studied in detail by Erdos and Renyi (E-R) 
m. For RGs, there is a hxed probability p for two nodes being connected 
and one can show that for a sufficiently large number of nodes N, the degree 
distribution P{k) tends to a Poissonian. 

The E-R model guided our thinking about complex networks for many decades 
until the discovery by Barabasi and co-workers uma nearly two decades 
back that the topology and structure of most networks around us are rad¬ 
ically different. For example, many networks in the real world such as, the 
World Wide Web (WWW) [13], networks of social interactions [Tl|, protein 
and metabolic networks [15] and technological networks [16] , tend to obey cer¬ 
tain self-organizing principles in their evolution either due to some inherent 
property of the system or due to the nature of human interactions as in social 
networks. The topology of such networks shows a scale invariance with the 
degree distribution obeying a power law, P{k) oc and the value of 7 is 
found to be typically between 2 and 3. The discovery of such scale-free (SF) 
networks triggered a lot of interest in the theory of complex networks. 

Along with the discovery of SF networks, the concept of small-world (SW) 
networks was introduced by Watts and Strogatz Ha- Though the classical RGs 
are amenable to a great deal of mathematical analysis, they are poor models 
as far as real networks are concerned. Firstly, they show poor clustering and 
their clustering coefficient (GG) is directly proportional to p. Secondly, for a 
given p, as the number of nodes N increases, the average degree < k > also 
increases correspondingly. Gonsequently, for p above a threshold value, the 
characteristic path length (GPL) of RGs remains very small and independent 
of N. 

Watts and Strogatz (W-S) showed that, starting from a ring lattice of N nodes 
with nearest neighbour coupling and randomly re-wiring a small fraction fd of 
edges, results in a complex network with high GG compared to the RG and 
small GPL comparable to that of a RG for a range of values of (d. Moreover, 
for such networks, as N increases, the GPL increases only as log N. Thus 
the W-S model displays many characteristic properties of real world networks 
and provides one possibility of obtaining the SW property, often found in real 
world networks, with different levels of complexity by tuning fd. 

The above developments resulted in complex networks and the related mea- 


2 




sures being applied as tools in many areas of research. The most recent among 
them is the analysis of time series from dynamical systems using statistical 
measures of complex networks. To study many dynamical processes in the real 
world, one often resorts to the analysis of time series obtained from the system. 
An area of special interest is where the underlying system appears to show de¬ 
terministic nonlinear behavior. The methods and measures of nonlinear time 
series analysis and chaos theory [T8IIT9] are commonly used in such cases. In 
the last few years, measures based on complex networks have gained a lot 
of importance in nonlinear time series analysis. All such measures propose a 
mapping from the time series domain to the network domain and then proceed 
to characterize the dynamical system in terms of the statistical measures of 
the resulting complex network. By doing this, one expects to resolve compli¬ 
mentary features that are not captured by conventional methods of time series 
analysis, especially the structural and topological properties of the underlying 
chaotic attractor. 

Even though several methods have been proposed [^211122] to convert time 
series into networks, an approach incorporating the generic property of re¬ 
currence [23] of a dynamical system has been prominently applied for the 
conversion of time series into networks. In this method, the time series is hrst 
embedded in a suitable dimension M using time delay co-ordinates [23] to 
reconstruct the attractor. Every point on the attractor is identihed as a node 
and the network can be constructed in two ways, either by taking a hxed 
number of nearest neighbours [2S] or by taking a hxed hyper-sphere of radius 
e with the point as the centre. In this work, we consider the second method 
for the construction of the network where a reference node i is connected to 
another node j if the Euclidean distance dij between the corresponding points 
on the attractor in the reconstructed space is less than or equal to the recur¬ 
rence threshold e, that is, if dij < e. The resulting complex network, called the 
e - recurrence network or simply recurrence network (RN) [231127] . has been 
shown to have great potential for a wide range of practical applications, from 
identifying critical transitions in dynamical systems [28] to the classihcation 
of cardio-vascular time series [29]. Note that, by construction, the RN is an 
undirected and unweighted graph with a symmetric and binary adjacency ma¬ 
trix A, with elements = 1 or 0 , depending on whether the two nodes i and 
] are connected or not. 

Though RNs and related statistical measures are widely applied in nonlinear 
time series analysis, their properties, as the threshold is increased, are not 
fully understood from a complex network perspective. It is well known that 
all RNs have two properties in common. Firstly, the degree distribution of 
every RN is unique and is closely related to the probability density variations 
over the embedded attractor from which it is mapped [30] . We will discuss this 
in detail below. Secondly, there is an absence of long range connections as the 
maximum edge length is limited by the recurrence threshold e. By dehnition. 
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Fig. 1. Construction of the RN from the time series of the standard Rossler at¬ 
tractor’s y - component using time step At = 0.05 with time delay r = 24 and 
M = 3. The time series and the embedded attractor are shown (top and bottom, 
respectively) in the left panel while the RN and its degree distribution are shown 
in the right panel. The RN is constructed taking every point on the attractor as a 
node and connecting every node to all other nodes within a recurrence threshold 
of e = 0.1 (see text). Error bars of P{k) originate from the fact that for a network 
with N nodes, the number n of nodes with a given degree k has a standard error of 
y/n{k). For n{k) —)■ 0, its value is normalised to 1, the minimum count. 


RNs are random geometric graphs (RGG) in the considered system’s phase 
space [3n|31ll32j . Here RGGs are RGs where each vertex is randomly assigned 
co-ordinates in some metric space according to some prescribed probability 
distribntion fnnction, and vertices are connected if and only if they are sepa¬ 
rated by less than a certain maximnm distance |33j . 

In this Letter, we nnmerically investigate the specihc properties of RNs using 
three primary measures of a complex network, the degree distribution, the GG 
and the GPL. In particular, we consider the range of threshold values beyond 
the small operational window usually used for the construction of the RN and 
check whether the resulting network can show the properties of either RGs or 
SF networks or behave as a small-world network. 
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2 Numerical Results 


All numerical simulations are done using N = 2000 nodes, unless otherwise 
mentioned. For the construction of SF networks, we use the specihc scheme of 
preferential attachment discussed in detail in [M] . Here we start with a small 
number of nodes mo- A new node is then added which gets connected to m 
number of existing nodes. This process is repeated, increasing the number of 
edges by m for each newly connected node. By changing either mo or m we 
can construct SF networks with different 7 . Here we do both and construct 
the SF networks using 3 values for mo (2, 4 and 10) and in each case, use 3 
different values for m, namely, 1, 2 and 4. Moreover, the RGs are constructed 
for different values of p. Time series from several standard low-dimensional 
chaotic systems are used for the construction of RNs. For continuous systems 
(flows) in 3D, we use the embedding dimension M = 3 and for discrete systems 
(maps) in 2D, we use M = 2. For all systems, we use the time series from the 
y - component and for all continuous systems, the time step At used for 
numerical integration is 0.05 with the time delay computed by the automated 
algorithmic scheme proposed by us [55] . The values of time delay used for the 
Rossler, Lorenz, Duffing and Ueda systems are 24, 6 , 25 and 13 respectively. 

To get a quantitatively comparable value for the percolation threshold e and 
to make the comparison between systems possible, we hrst transform the time 
series to a uniform deviate so that the size of the reconstructed attractor is 
re-scaled into a unit cube in M dimensions. This transformation is not a trivial 
rescaling as it stretches the embedded attractor uniformly in all directions. We 
have already shown [351136] how effective this transformation is in computing 
conventional measures such as, correlation dimension D 2 and entropy K 2 , for 
low as well as high-dimensional systems from time series [57]. As a result 
of the uniform deviate transformation, we have found that it is possible to 
have a small identical range of recurrence threshold Ae, that can be taken as 
the operational window for constructing the RN from time series of different 
systems for a given embedding dimension M [38]. Here we choose the value 
of e as the minimum of this range where the giant component of the RN just 
appears, as suggested by Donges et al. [52]. The value of e is found to be 
0.06 for M = 2 and 0.1 for M = 3 for N < 10000 as discussed in detail 
in [38]. The RN from a random (white noise) time series is also constructed 
using M = 3 (for comparison with the RNs from chaotic systems) whose 
degree distribution is Poissonian with < k >~ 7 for the selected e. We hnd 
that this distribution coincides with the degree distributions of a RG for 2000 
nodes with p = 0.0035 and a RGG in 3 dimensions with the spatial range 
of connections limited to the same threshold e = 0.1 as that of the RN. To 
draw the connection with the percolation threshold of the RN, one can refer 
to [52]. In Fig. [H we show the construction of a typical RN from a standard 
Rossler attractor time series. The node degrees of the RN vary over a wide 
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Fig. 2. Cumulative degree distribution of RNs from several low-dimensional chaotic 
attractors and a random time series. The chaotic attractors are the Henon (open 
square), Lozi (solid square), Lorenz (open circle), Rossler (solid circle). Duffing 
(open triangle) and Ueda (solid triangle) systems whose D 2 ranges are also indi¬ 
cated. The parameters used for all the systems are given in [39]. We observe that 
the position where F{k) decreases sharply follows a pattern in accordance with the 
dimension of the attractor. 

range determined by the local probability density variations over the attractor. 
Specihcally, the topology of the RN closely resembles that of the embedded 
attractor, with nodes corresponding to regions of high probability density in 
the attractor having higher degree in the RN and vice versa. This is evident 
from Fig. [H 
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Fig. 3. CPL-CC combined plot for different classes of networks with 2000 nodes in 
all cases except the W-S network. We show the results for 3 SF networks (solid 
triangles) with 7 = 2.28,2.49,2.67, 3 RGs (solid circles) with link probability 
p = 0.002,0.0035,0.02 and 3 W-S small-world networks (solid squares) with 100 
nodes (to reduce computation time and memory) using /3 = 0.04,0.05,0.06. In all 
cases, the size of the symbols increases with the values of 7 , p and /3 respectively. 
The open star on the right side corresponds to the RN from a random time series 
and the asterisk corresponds to a RGG with the range of connection limited to that 
of the RNs in 3 dimension, namely, e = 0.1. These two networks are completely 
different from the other networks in terms of GPL-CC values. 
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As the dynamical system evolves, the probability density over the attractor is 
preserved. The degree fcj of a reference node i represents the local connectivity 
of the RN and corresponds to the local phase space density aronnd the refer¬ 
ence point in the attractor. It has already been shown [2^11^ that the local 
probability density variation is reflected in the degree distribntion of the RN: 


N 


oc p(fi 


( 1 ) 


for a given M, where p(?y) is the invariant density aronnd the point r*. In the 
light of this eqnation, ki/N can be called the degree density [31]. The above 
eqnation tells ns that the distribntion in terms of the degree density remains 
invariant for a given attractor as has been explicitly shown by us earlier [38] . 

Due to this close connection with the invariant density, the degree distribution 
also carries important structural information on the corresponding chaotic at¬ 
tractor in a subtle manner. To show this explicitly, we consider the cumulative 
degree distribution of the RNs from several low-dimensional chaotic attrac¬ 
tors. This representation is often used in the analysis of power-law degree 
distributions to retrieve the scaling that is not visible in the power-law tail. 
It has been extensively discussed in the literature [2] in connection with the 
analysis of SF networks and is dehned as 

^max 

F{k) = P(k') (2) 

k'=k 


It can be shown that if P{k) is a power-law with exponent 7 , then F{k) will 
also follow a power-law with exponent 7 — 1. Here we make use of F{k) only 
to get a better comparison between the degree distributions of various chaotic 
attractors. 

In Fig. 121 we show the variation of log F{k) as a function of log k for RN 
from several low-dimensional chaotic attractors and random time series. The 
parameters used for the generation of all the chaotic attractors are given in the 
Appendix of [3H]- In all cases, we use the y - component of the system with 
time step 0.05 for continuous systems. Note that the position of the tail of 
the distribution follows a pattern in accordance with the dimension D 2 of the 
attractor. To get a possible explanation for this observation, we have studied 
the actual distribution F{k) of the RNs. In all cases, there is a range of k 
values and a particular k value, say k*, at which F{k) has a maximum value. 
The cumulative distribution starts decreasing sharply as k exeeds k*. Now, 
the value of k* increases with the range of k values which, in turn, depends on 
the clustering of points on the attractor, as the attractor is converted into the 
network. For Henon and Lozi attractors, the clustering is maximum as they 



are characterized by a lower D 2 value and consequently the range of k values 
is maximum. This makes their prohle getting shifted to extreme right. As D 2 
increases, the range of k values and hence k* decreases due to the rescaling 
of all systems to the unit cube, shifting the tail to lower k values. For the 
random time series, when computed from a hnite time series with embedding 
dimension M, we get D 2 = M and the range is minimum with k* =< k > and 
hence the tail of F{k) is located at the extreme left. The above argument is 
supported by numerical results from earlier studies [26] , if the local clustering 
coefficient can be understood analytically from the theory of random geometric 
graphs as being associated to a “mean local clustering dimension” dehned in 
|31] . which should behave similar to the classical correlation dimension D 2 
considered here. 


Next, we consider how the RNs are different with respect to the other two 
measures, CC and CPL. We generate an ensemble of SF networks with differ¬ 
ent 7 values using the specihc scheme of preferential attachment as mentioned 
above and RGs with different link probability p. The CPL is dehned by the 
equation 


< I >= 


1 

N{N -1) 






( 3 ) 


where kj is the shortest path length for all pairs of nodes {t, j) in the network. 
The CC of the network is dehned through a local clustering coefficient c„ 
which measures the probability that two randomly chosen neighbours of a 
given vertex v are mutually linked. For hnite graphs, this is estimated [31] in 
terms of relative frequency of links between the neighbouring vertices for the 
node V. 


C71 — 
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The average value of Cy is taken as the CC of the whole network: 


CC 



Cy 


( 5 ) 


The results of our computations are shown in Fig. [3] using a combined CPL- 
CC graph for N = 2000 in all cases except the W-S network. We have shown 
the values for the RN of a random time series and those from 3 diherent SF 
networks, RCs and W-S small-world networks. The values for a RCC in 3 
dimensions with the range limited to that of the RN (in 3 dimensions) are 
also shown. Two results are evident from the hgure. The CPL values of all 
networks except the RN and RCC are very small and always < 6 . Secondly, 
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Fig. 4. CPL-CC values for RNs generated from several low-dimensional chaotic at¬ 
tractors and random time series for 4 different values of N, namely, 1000, 2000,5000 
and 10000 with the point denoted as 1 representing N = 1000. For comparison, 
values for a RG (with p = 0.004) and a SF network (with 7 = 2.162) with the same 
N values are also shown. 

the RG with p = 0.0035, whose degree distribution coincides with that of the 
RN from the random time series behaves completely differently in the CPL-CC 
graph. The reason for this discrepancy is the absence of long range connections 
in the RN. In fact, this is true for all RNs and this property cannot change 
with the increase in the number of nodes N. Both CPL and CC are found 
to converge and saturate to a finite value for all RNs and cannot fall below a 
certain limiting value which depends on e. This limiting value is much higher 
compared to the typical values for SF networks and RGs. This is shown in 
Fig. musing RNs constructed from several standard chaotic attractors with N 
varying from 1000 to 10000. Our numerical results show that RNs, constructed 
with e within the usual operational window, cannot show SW property. 

Suppose that we now increase the value of e above the recurrence threshold 
which, in turn, increases the range of spatial connection for each node in 
the RN. The resulting RN may not be able to characterize the statistical 
properties of the embedded attractor, but one can look at the properties of 
the resulting network from a complex network perspective. For example, in the 
case of spatial networks, such as, communication and transportation networks, 
Barthelemy [iOlITT] has shown that as the interaction range between nodes is of 
the order of system size or larger, the spatial effects become negligible and the 
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Fig. 5. Variation of the CPL and CC of RNs when the recurrence threshold e is 
increased from a value close to the percolation threshold e = 0.1 with M = 3. The 
RN first undergoes a smooth transition to a network with high CC (~ 0.8) and 
small CPL (< 5), capable of showing SW property for an intermediate range of e 
approximately between 0.25 and 0.5. As e —>■ 1, the RN is found to cross over to 
the classical RG (in the limit p —>■ 1), where both CC and CPL —)• 1 independent 
of N. We also show the variation of CC and CPL for a W-S network as a function 
of the re-wiring probability /3 (solid line) and the same for a RG with connection 
probability p (dashed line). The number of nodes in all cases is 2000. 
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Fig. 6. Variation of CPL as a function of log N for RNs corresponding to the 3 
ranges of threshold e used for the construction of the network. The top panel (a) 
corresponds to e within the range of recurrence thresholds where the RN represents 
the characteristic properties of the attractor and the CPL converges to a limiting 
value as N increases. The middle panel (b) corresponds to intermediate values of e 
(between 0.25 and 0.5) where the RN is capable of displaying the SW property as 
CPL increases linearly with log N. The bottom panel (c) corresponds to large values 
of e (> 0.5) where the RN has long range connections and the CPL tends to a small 
value independent of N. For (b), the value of e is tuned within the respective range 
keeping < k > constant as N is increased. We find the SW property for a range of 
< k > values between 600 to 1500 depending on the value of e. The graph in (b) is 
plotted with e = 0.3 and decreasing e as is increased by fixing < k >= 800. In all 
three cases shown above, the solid triangles connected by a dashed line represent the 
CPL of Lorenz attractor while the solid circles connected by a solid line represent 
that from random time series. 
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networks tend to show SF property. RNs are spatially constrained networks 
and hence we now check whether, under some limiting conditions, the RNs 
can show properties of either SW networks or classical RGs. 


In Fig. m we show the variation of CC and CPL for the RNs from random time 
series and Lorenz attractor as a function of e. While the CC shows a steady 
increase with e, the CPL decreases much more quickly initially and saturates 
to the minimum possible value of 1. We also show the CC and CPL for the 
RC (dashed line) for increasing values of p and a typical Watts-Strogatz SW 
network (solid line) as a function of increasing re-wiring probability f3. We hnd 
that for an intermediate range of e values, approximately from 0.25 to 0.5, the 
RNs show high CC (~ 0.8) and small CPL (< 5) which are the basic criteria 
for a network to show the small-world property. To check this, we increase 
the value of N keeping the average degree < k > approximately constant by 
changing the value of e within the above range ( between 0.25 and 0.5) and 
compute the value of CPL as function of log N. The result, presented in Fig. [B] 

(b) , clearly shows the small-world property. Thus, the intermediate range of e 
can be considered as a small-world phase for the RNs. As e further increases 
and tends to the system size, all RNs smoothly cross over to the classical 
RCs with high p where the CPL and CC —)■ 1 independent of N. To get a 
comparison of the behaviour between different regimes of e, we also show the 
variation of CPL with log N for RN with e equal to the recurrence threshold 
in Fig. E] (a) and for the case of long range connections with e > 0.5 in Fig. O 

(c) . 


We have checked this result for two more standard chaotic attractors, namely, 
the Rossler attractor and the Duffing attractor and found an identical behav¬ 
ior. Hence, based on our numerical results, we conjecture that all RNs pass 
through an intermediate small-world phase as the value of e is increased, before 
undergoing a transition to the classical RC for large e. Our result provides a 
method for constructing a SW network from time series of dynamical systems 
by adjusting the range of connections between the nodes. The result is also 
relevant to search for the SW property in network models that require a multi¬ 
dimensional system with a metric, as for example in the case of spreading of 
diseases [12], in continuous percolation theory [H] and in the study of random 
walks on fractals [Hj. For example, it is well known [33] that the spreading 
of diseases can be modeled by a RCC. If the nature of spreading is limited to 
nearest neighbour interaction (say, controlled by contact only), the network 
cannot show SW property. However, if the spreading occurs by other means, 
so that the range e is large, the network can very well switch over to a SW. 
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3 Conclusion 


In this Letter, we analyse the RNs, which are unweighted and undirected 
graphs, from a complex network perspective. Three primary characteristics of 
a network, namely, the degree distribution, the CC (an averaged local measure) 
and the CPL (a global measure) are used for the analysis. We hnd that the 
degree distribution of every RN is unique and characteristic to the probability 
density variations over the representative attractor. By computing the CC and 
the CPL of several RNs from standard chaotic attractors, we explicitly show 
that the RNs constructed with the recurrence threshold at and just above the 
percolation threshold are basically different with both CC and CPL being high 
due to the absence of long range connections between the nodes. 

The characteristic properties of the RN change with the value of e used to 
construct it. As e —)■ 0, CPL —)■ oo and CC —)■ 0. For values of e below the 
recurrence threshold, the RN consists of several disconnected clusters. For e 
within the operational window above the percolation threshold, the network 
becomes globally connected (with the appearance of a giant component) and 
the RN measures characterize the statistical properties of the underlying sys¬ 
tem. In this phase, CPL remains relatively high (with the network unable 
show SW property) and converges to an asymptotic value as N is increased. 

When the value of e is increased further, the constraints induced by the short 
connectivity range slowly disappear with the emergence of long range connec¬ 
tions and all the RNs initially cross over to a SW phase with high CC and 
small CPL that increases as log N. Finally, as the value of e ^ 1, the RN 
undergoes a smooth transition to classical RGs with p ~ 1 with both CPL 
and CC —)■ 1 independent of the number of nodes N. Thus, for all RNs, there 
are basically three phases in terms of the range of e values used for the con¬ 
struction of the RN, namely, a small operational window of e that captures 
the statistical properties of the corresponding attractor, an intermediate range 
where the RN can display the SW property and hnally e of the order of the 
size of the attractor where the RN crosses over to the classical RG. 
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